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Abstract. In this talk I will introduce the principle of stochastic stability and dis- 
cussing its consequences both at equilibrium and off-equilibrium. 



INTRODUCTION 

In this talk I will underline the physical meaning of replica symmetry breaking 
[1-3]. In this framework I will introduce the principle of Stochastic Stability: I will 
present its far reaching consequences on the equilibrium properties of the system 
and on the off-equilibrium behaviour. 



THE COEXISTENCE OF MANY PHASES 

Usually, if a system has different phases which are separated by a first order tran- 
sition, just at the phase transition point a very interesting phenomenon is present: 
phase coexistence. This usually happens if we tune one parameter. This behaviour 
is summarized by the Gibbs rule which states that, in absence of symmetries, we 
have to tune n parameters in order to have the coexistence of n + 1 phases. 

In the case of complex systems we have that the opposite situation is valid: the 
number of phases is very large for a generic choice of parameters. It is usual to 
assume that all these states are globally very similar: states can be separated only 
by comparing one state with an other. An example of this phenomenon would be 
a very long heteropolymer, e.g, a protein or RNA, which folds in many different 
structures. However quite different foldings may have a very similar density. Of 
course you will discover that two proteins have folded in two different structures if 
we compare them. 

In order to be precise we should consider a large but (finite) system [3] . We want 
to decompose the phase state in valleys (phases, states) separated by barriers. If 
the free energy as function of the configuration space has many minima, the number 
of states will be very large. 



Let us consider for definitiveness a spin system with N points (spins are labeled 
by i, which in some cases will be a lattice point). States (labeled by a) are char- 
acterized by different local magnetizations: m a (i) = (a(i)) a , where (-) a is the 
expectation value in the valley labeled by a. The average done with the Boltz- 
mann distribution is denoted as (■) and it can be written as linear combinations of 
the averages inside the valleys. We have the relation: 

(•>«E w «0«- (!) 

a 

We can also write the relation w a oc exp(— /3F a ), where by definition F a is the free 
energy of the valley labeled by a. 

In the rest of this talk I will call J the control parameters of the systems. The 
average over J will be denoted by a bar (e.g. F). I will consider here the case 
where a quenched disorder is present: the variables J parametrize the quenched 
disorder. 



THE OVERLAP AND ITS PROBABILITIES 

We have already remarked, states may be separated making a comparison among 
them. At this end it is convenient to consider their mutual overlap. Given two 
configurations, we define their overlap: 

?M = 1 £ a(i)r{i). (2) 

iV i=l,N 

The overlap among the states is defined as 

g(a,7) = T7 Z m^m-Xi) ~ q[<T,r] , (3) 

ly i=l,N 

where a and r are two generic configurations that belong to the states a and 7 
respectively. 

We define Pj{q) as probability distribution of the overlap q at given J, i.e. the 
histogram of q[a, r], where a and r are two equilibrium configurations. Using eq. 
(1), one finds that 

P-M) = ^w a w 1 5(q - q an ) , (4) 

where in a finite volume system the delta functions are smoothed. If there is more 
than one state, Pj(q) is not a single delta function: 

Pj(q) + S(q - q EA ) ■ (5) 

If this happens we say that the replica symmetry is broken: two identical replicas 
of the same system may state in a quite different state. However at zero magnetic 




FIGURE 1. The function P(q) after average over many samples in four dimensions for system 
of size i 4 , (L=3. . . 10) . 

field Pj(q) is an even function, so that at low temperature it will contains two delta 
functions at \q\ ^ 0. 

There are many models where the function Pj(q) is non trivial: a well known 
example is given by Ising spin glasses [1,4,5]. In this case the Hamiltonian is given 
by 

H = - J i,kVi<7k ~ h i a i ' ( 6 ) 

i,k i 

where a = ±1 are the spins. The variables J are random couplings (e.g. Gaussian 
or ±1) and the variables hi are the magnetic fields, which may be point dependent. 
Let is consider two different models for spin glasses: 

a) The Sherrington Kirkpatrick model (infinite range): all N points are con- 
nected: Jj ; fc = 0(N^ 1 ^ 2 ). Eventually N goes to infinity. 

b) Short range models: % belongs to a L D lattice. The interaction is nearest 
neighbour (the variables J are or zero or of order 1) and eventually L goes to 
infinity at fixed D (e.g. D = 3). 

Analytic studies have been done in the case of the SK model, where one can prove 
rigorously that the function Pj(q) is non-trivial. In the finite dimensional case a 
homologous theorem has not been proved and in order to answer to the question if 
the function Pj(q) is trivial we must resort to some numerical simulations [6] or to 
experiments. 

In this case numerical simulations present ample evidence that in three and four 
dimensions the function Pj(q) is non-trivial, i.e. at zero magnetic field it is not the 
sum of two delta functions: this can be seen looking to the average of Pj(q), shown 
fig. (1) from [7]. We will see later that the the function Pj(q) strongly depends on 
J. 
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FIGURE 2. The quantity (g 2 ) 2 - (±(g 4 ) + |(g 2 ) ) as function of the temperature for different 
values of L in D = 3. 

STOCHASTIC STABILITY 

In the nutshell stochastic stability states that the system we are considering 
behaves like a generic random system [8-10]. Technically speaking in order to 
formulate stochastic stability we have to consider the statistical properties of the 
system with Hamiltonian given by the original Hamiltonian (if) plus a random 
perturbation (Hr): 

H(e) = H + eH R . (7) 

Stochastic stability states that all the properties of the system are smooth functions 
of e around e = 0, after the appropriate averages over the original Hamiltonian and 
the random Hamiltonian. 

Typical examples of random perturbations are (we can chose the value of r in an 
arbitrary way): 

R (r) = N (r-l)/2 £ R{n . . . ir)(j{il) . . . a{ir) ^ (8) 
i\...i r 

where for simplicity we can restrict ourselves to the case where the variables R are 
random uncorrelated Gaussian variables. 

It useful to remark that if a symmetry is present, a system cannot be stochasti- 
cally stable. Indeed spin glasses may be stochastically stable only in the presence of 
a finite, non zero magnetic field which breaks the a <-> — a symmetry. If a symme- 
try is present, stochastic stability may be valid only for those quantities which are 
invariant under the action of the symmetry group. It is also remarkable that the 
union of two non-trivial uncoupled stochastically stable systems is not stochastically 



stable. Therefore a non-trivial stochastically stable system cannot be decomposed 
as the union of two or more parts whose interaction can be neglected. 
In the general case stochastic stability implies that 

P(qi, to) = PMWM) = \p{Qi)P{Q2) + \p(qi)5(qi - q 2 ) • (9) 
A particular case of the previous relation is the following one: 

'W = l~¥) + l'¥) 2 ■ (io) 

We have tested the previous relations in three dimensions as function of the 
temperature at different values of L [3]. In fig. 2 we plot the quantity (q 2 ) 2 — 

(f (o' 4 ) + \{q 2 ) )) which should be equal to zero. Indeed it is very small and its 
values decreases with L. The two quantities in eq. (10) in the low temperature 
region are a factor of 10 3 bigger of their difference. I believe that there should be 
few doubts on the fact that stochastic stability is satisfied for three dimensional 
spin glasses. 

OFF-EQUILIBRIUM DYNAMICS 

The general problem that we face is to find what happens if the system is carried 
in a slightly off equilibrium situation. There are two ways in which this can be 
done. 

a) We rapidly cool the system starting from a random (high temperature) con- 
figuration at time zero and we wait a time much longer than the microscopical one. 
The system orders at distances smaller that a coherence distance £ (t) (which even- 
tually diverges when t goes to infinity) but remains always disordered at distances 
larger than 

b) A second possibility consists in forcing the system in an off-equilibrium state 
by gently shaking it. This can be done for example by adding a small time de- 
pendent magnetic field, which should however strong enough force a large scale 
rearrangement of the system [11]. 

In the first case we have the phenomenon of ageing. This effects may be eviden- 
ziated if we define a two time correlation function and two time relaxation function 
(we cool the system at time 0) [12,13]. The correlation function is defined to be 

1 N 

c(t, t w ) = — J2(^(t w )^(t w + 1)) , (ii) 
iV i=i 

which is equal to the overlap q(t w , t w + t) among a configuration at time t w and one 
at time t w + t. The relaxation function S(t, t w ) is a just given by 



(12) 
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FIGURE 3. Relaxation function versus correlation in the Edwards- Anderson (EA) model in 
D = 3 T = 0.7 ~ jT c and theoretical predictions from eq. (16) . 



where 8m is the variation of the magnetization when we add a magnetic field 
5h starting from time t w . More generally we can introduce the time dependent 
Hamiltonian: 



H = H + 6{t-t w )J2hi<Ji 



(13) 



The relaxation function is thus defined as: 



We can distinguish two situations: 



1=1 



dhi 



(14) 



• For t « t w we stay in the quasi-equilibrium regime [14], C(t,t w ) ~ C eq (t), 
where C eq (t) is the equilibrium correlation function; in this case qEA = 
lim^oo lim^^oo C(t, t w ). 

• For t = 0(t w ) or larger we stay in the aging regime. If simple aging holds 
C(t,t w )(xC(t/t w ). 

In the equilibrium regime, if we plot parametrically the relaxation function as 
function of the correlation, we find that 



dS_ 

dC 



-1 



(15) 



which is a compact way of writing the fluctuation-dissipation theorem. 



Generally speaking the fluctuation-dissipation theorem is not valid in the off- 
equilibrium regime. In this case one can use stochastic stability to derive a re- 
lation among statics properties and the form of the function S(C) measured in 
off-equilibrium experiments [12,13,15]: 

-% = £dqP(q) = X(C). (16) 

The validity of these relation has been intensively checked in numerical experi- 
ments (see for example fig. (3) from [16]). 

In spin glasses the relaxation function has been experimentally measured many 
times in the aging regime, while the correlation function has not been measured: 
it is a much more difficult experiment in which one has to measure thermal fluc- 
tuations. Fortunately enough measurements of both quantities for spin glasses are 
in now progress. It would be extremely interesting to see if they agree with the 
theoretical predictions. 
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